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ABSTRACT 

The cosmic web that characterizes the large-scale structure of the Universe can 
be quantified by a variety of methods. For example, large redshift surveys can 
be used in combination with point process algorithms to extract long curvi-linear 
filaments in the galaxy distribution. Alternatively, given a full 3D reconstruction of 
the velocity field, kinematic techniques can be used to decompose the web into voids, 
sheets, filaments and knots. In this paper we look at how two such algorithms - the 
Bisous model and the velocity shear web - compare with each other in the local 
Universe (within 100 Mpc), finding good agreement. This is both remarkable and 
comforting, given that the two methods are radically different in ideology and applied 
to completely independent and different data sets. Unsurprisingly, the methods are in 
better agreement when applied to unbiased and complete data sets, like cosmological 
simulations, than when applied to observational samples. We conclude that more 
observational data is needed to improve on these methods, but that both methods 
are most likely properly tracing the underlying distribution of matter in the Universe. 
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1 INTRODUCTION 

That galaxies are not homogeneously and uniformly dis¬ 
tributed in the Universe has been known at least since 
Reynolds (1923) stated that "... there is very definite ev¬ 
idence for a band [of galaxies], fairly widespread ...” (see 
also Herschel 1847), a feature later confirmed and named 
the “supergalactic plane”. In modern terms, the anisotropic 
collapse of density perturbations into large plane-like pan¬ 
cakes was predicted in the seminal work of Zel’dovich (1970). 

The concepts of filaments in the galaxy distribution was first 
suggested by Joeveer et al. (1978), while Bond et al. (1996) 
argued that the large-scale distribution of matter constitutes 
a network known as the cosmic web. 

Since then, many researchers have sought to quantify 


cosmology: theory - dark matter - large-scale 

the cosmic web using a variety of techniques to gain in¬ 
sight into the relationship between galaxies and their en¬ 
vironment. Often simulations are used as a testbed, since 
the full 6 dimensional phase space of the matter distribu¬ 
tion is available. Sophisticated algorithms that attempt to 
capture the cosmic web’s multi-scale nature (e.g. Cautun 
et al. 2013), or that are based on shell crossing and phase- 
space folding (Falck et al. 2012; Shandarin et al. 2012; Abel 
et al. 2012) have proven useful to understand how the cos¬ 
mic web evolves (see also Aragon-Calvo et al. 2007, 2010). 
Other methods, for example based on the Hessian of the tidal 
(Hahn et al. 2007; Forero-Romero et al. 2009) or shear ten¬ 
sor (Hoffman et al. 2012; Libeskind et al. 2014) have proven 
useful to study aspects of haloes such as their mass function 
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(Metuki et al. 2015) or spin (Libeskind et al. 2012; Codis 
et al. 2012). 

The construction of the cosmic web from astronomical 
data represent a considerable challenge. A partial list of ob¬ 
stacles to overcome includes magnitude limits, volume lim¬ 
its, biases of multiple sorts, incompleteness, redshift-space 
distortions, obscuration by the Galactic disc and sparse sam¬ 
pling. Nevertheless, various techniques to find filaments in 
observational samples have been devised (e.g. Sousbie et al. 
2008, 2011; Bond et al. 2010; Stoica et al. 2010; Alpaslan 
et al. 2014; Hong & Dey 2015, among others). The cor¬ 
respondence between numerical and observational methods 
remains unclear and further study is clearly necessary. 

As an attempt to resolve this ambiguity, Tempel et al. 
(2014) examined the alignment between cosmic filaments de¬ 
tected by a point process with those detected by exploiting 
the velocity shear tensor, using an A-body simulation. They 
found that indeed these methods agree to a large extent. 
The work of Tempel et al. (2014) is extended here from 
simulations and applied to two data sets representing the 
large-scale structure of the local (<100 Mpc) Universe. We 
compare the filaments constructed from the Bisous method 
applied to the 2MRS (Huchra et al. 2012; Lavaux & Hudson 
2011) catalogue of galaxies, to those obtained by the velocity 
shear method (Hoffman et al. 2012) from the Wiener Filter 
reconstruction of the Cosmicflows-2 (CF2, Tully et al. 2013) 
data base. 

The agreement that Tempel et al. (2014) found in sim¬ 
ulations indicated that these two methods trace the same 
cosmic structures. However a comparison between these two 
methods applied to simulations is bound to give a better 
agreement than when applied to observational data. For 
one, Tempel et al. (2014) had the full density and velocity 
field at their disposal from which to compute the shear field 
and Bisous filaments; when using observations only magni¬ 
tude limited, partially obstructed, galaxy redshift surveys 
are available for the latter and sparse, incomplete and in¬ 
homogeneous catalogues of radial velocities are available for 
the former. In Tempel et al. (2014), all dark matter haloes 
(above a minimum mass cut) in the computational box were 
used as proxies of the galaxy distribution returning fairly 
complete filaments throughout the simulated volume; in ob¬ 
servations, galaxies obtained according to the survey’s se¬ 
lection function is used for this purpose. Since simulations 
offer the full 6D phase space information of the density and 
velocity fields to high accuracy, the shear tensor can be con¬ 
structed on very small scales. When dealing with the sparse 
sampling of the 4D (position and radial velocity) phase space 
information of observations, the resolution of the shear ten¬ 
sor computation is significantly worse. It is unclear whether 
these two methods can agree at all on observational data, 
given the long list of limitations. 


2 METHOD 

2.1 V-web filaments from Cosmicflows-2 

The CF2 data (Tully et al. 2013) are a comprehensive com¬ 
pilation of the universe’s radial peculiar velocity field, within 
around 100 Mpc of the Milky Way. The radial peculiar ve¬ 
locity field can be used to reconstruct the full 3D velocity 
and density fields by means of a Wiener Filtering algorithm 
(Courtois et al. 2012). We refer the interested reader to 
the important papers describing Wiener Filtering and con¬ 
strained realizations for details (Hoffman & Ribak 1991; La- 
hav et al. 1994; Zaroubi et al. 1995, 1999; Fisher et al. 1995) 
and explain the method only schematically here. 

The Wiener Filter is essentially a linear minimum vari¬ 
ance solution to recover an underlying, continuous field given 
a set of noisy, sparse, and incomplete data, such as the CF2 
survey. The technique assumes a prior for the form of the 
covariance matrix of the desired field. In this case, the prior 
is derived from the assumption of a ACDM universe where 
structures form out of Gaussian random fields. The Wiener 
Filtering algorithm takes as an input the measured radial ve¬ 
locity field and returns a reconstructed 3D (linear) velocity 
field which can be used to quantify the cosmic web kinemat¬ 
ically. There is an implicit smoothing in the Wiener Filter 
which is adaptive and inversely proportional to the quality 
of the data, namely the sampling and magnitude of the ob¬ 
servational errors. Where the observational errors are large 
or where the sampling is poor, the Wiener Filter has a large 
effective smoothing. In practice, the reconstructed velocity 
and density fields are attenuated with the depth, since ob¬ 
servational errors are large and sampling is poor at great 
distances. In the limit of vanishing data the Wiener Filter 
converges to the null fields. Nearby the effective resolution 
is ~ 5 Mpc. 

A kinematic classification can be made by examining 
the (reconstructed) velocity field, as described in Hoffman 
et al. (2012). A regular (256^) grid is placed on the velocity 
field and the shear tensor (similar in form to the “stress- 
strain” tensor used in mechanics) is computed and diago¬ 
nalized to return the principal components, known as the 
eigenvectors and corresponding eigenvalues. These describe 
the directions along which matter is principally collapsing 
(or expanding). The number of eigenvalues above an (arbi¬ 
trary) threshold can be used to classify a region as belonging 
to a void, sheet, filament, or knot. Typically the value for 
the threshold is chosen to recover the visual impression of 
the cosmic web (e.g. see Forero-Romero et al. 2009), but 
can also often taken to be null such that positive eigenval¬ 
ues correspond to axes along which compression occurs while 
negative eigenvalues imply expansion. In this work, we chose 
the latter threshold value of zero. 

With such a threshold, regions with 3 positive eigen¬ 
values (-I-, -b, -b) correspond to knots, 2 positive eigenval¬ 
ues (-b, -b, —) to filament, 1 positive eigenvalue (-b, —, —) to 
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sheets and 0 positive eigenvalues (—, —, —) to voids. Eigen¬ 
values and corresponding eigenvectors are conventionally de¬ 
noted as (Ai, Gi) with i — 1,2,3, respectively. In the case 
of filaments, one has a single negative eigenvalue (A 3 ); the 
corresponding eigenvector (es), defines the “spine” of the fil¬ 
ament. Its direction can be measured against the direction 
of filaments computed with the Bisous method, described 
below. 


2.2 Bisous filaments from marked point process 

Filaments in the local Universe can be found by linking 
galaxies from the 2MRS survey (Huchra et al. 2012; Lavaux 
& Hudson 2011) together in a curvi-linear fashion. The so- 
called “Bisous model” applied a marked point process with 
interactions (Stoica et al. 2005) to the 3D spatial distribu¬ 
tion of galaxies. The method provides a quantitative classi¬ 
fication that agrees well with a visual impression of the fil¬ 
amentary nature of the cosmic web. More details regarding 
the Bisous model can be found in Stoica et al. (2007, 2010) 
and Tempel et al. (2014). A brief, albeit simple summary is 
provided below. 

The Bisous model approximates the filamentary web by 
randomly placing small cylindrical segments on the galaxy 
distribution, and finding those cylinders which best “fit” the 
local, linear distribution of galaxies. The model assumes 
that locally the galaxy distribution can be approximately 
described by a series of relatively small cylinders, which can 
be combined to trace a filament if the neighbouring cylin¬ 
ders are oriented similarly. One of the advantages of such 
“point processes” is that they rely directly on the observed 
galaxy positions, and does not require any additional post¬ 
processing of the data for example the computation of con¬ 
tinuous density field via smoothing kernels. 

The solution provided by the Bisous model is stochas¬ 
tic and thus relies on the generation of Markov-chain Monte 
Carlo (MCMC). It is thus expected that there is some 
variation in the detected filamentary patterns for different 
MCMC runs of the model. However, thanks to the stochas¬ 
tic nature of the Bisous model, we gain a morphological and 
a statistical characterization of the filamentary pattern si¬ 
multaneously. 

The only free parameter in the model is the approximate 
radius of the filamentary cylinder. Once this is specific, the 
algorithm returns the filament detection probability field as 
well as the filament orientation field. Based on these two 
fields, a filament catalogue can be built in which every fila¬ 
ment is represented a specified axis, known as its spine. 

The spine detection follows two ideas. First, filament 
spines are required to be located where the algorithm speci¬ 
fies the filament detection probability is highest. In practice 
this is done by algorithmic examination of probability maps. 
Second, in these regions of high probability for the filamen¬ 
tary network, the spines are oriented along the orientation 


field of the filamentary network. See Tempel et al. (2014) for 
more details of the procedure. 

The filament scale in this study is chosen to be 
~ 0.7 Mpc. This scale and all other parameters in Bisous 
model are exactly the same as in Tempel et al. (2014). 
The chosen filament scale corresponds to the scale of galaxy 
groups/clusters and the detected filaments are the bridges 
between them. 


3 RESULTS 

After computing the Bisous filament catalogue from 2MRS 
and the direction (ea) corresponding to the spine of filaments 
extracted from the Wiener Filter reconstruction of the ve¬ 
locity field, we are able to measure the alignment of the two. 
Since the Wiener Filter and reconstructed velocity field per¬ 
forms best where the sampling is greatest (and attenuates 
to the mean field when the signal is lost), we restrict our 
analysis to regions of the universe that are within 100 Mpc. 
The fraction of the volume (within 100 Mpc) ascribed to 
knots, filaments, sheets and voids within this distance is 9.5, 
37.5, 41.9, 11.1%, respectively (shown also in Fig. 2, bottom 
panel). 

In practice, we find the locations of filaments in the 
Bisous catalogue of 2MRS galaxies first. Then, at the loca¬ 
tion of each filament, we select the corresponding ea eigen¬ 
vector from the reconstructed velocity field. The two are dot¬ 
ted and the distribution of the (cosine) of the angle between 
them is shown in Fig. 1. The two vectors are well aligned for 
all regions in the universe that can be probed (i.e. are as¬ 
sociated to Bisous filaments), as shown by the green line in 
each panel in Fig. 1. This indicates that Bisous filaments ap¬ 
pear to trace the cosmic velocity field well (and vice versa). 
Note that the two directions reported here originate from 
two vastly different data sets with different inherent biases: 
the Bisous filaments from galaxy redshift surveys and the 
V-web from reconstructions of the over-density and velocity 
field based on an opportunistic and incomplete compilation 
of galaxy distances using a variety of methods. Their close 
alignment indicates that these two completely independent 
methods are properly tracing the underlying distribution of 
matter. 

Although by definition, all locations that can be tested 
for alignment are in Bisous filaments, these regions of the 
Universe are not necessarily defined as filaments by the 
shear analysis. Indeed, these locations have a range of den¬ 
sities and kinematic classifications. In total, 16.9% of Bisous 
filament locations are V-web knots, 60.1% are V-web fila¬ 
ments, 21.1% are V-web sheets and 1.9% are V-web voids 
(shown also in Fig. 2, top panel). These fractions can be con¬ 
trasted with those found in the full volume to get a feeling 
for the nature of V-web defined environments that are sub¬ 
selected from the full reconstructed volume by the Bisous 
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cos(e) 

Figure 1. The probability distribution of the cosine of the angle 
formed between Bisous filaments and 03 , the axis of slowest col¬ 
lapse or fastest expansion computed from the velocity shear ten¬ 
sor. In each panel, the green dashed line indicates the alignment 
for all Bisous filaments. The Bisous filaments are then subdivided 
by V-web classification into voids, sheets, filaments, and knots 
(top to bottom) and the probability distribution in each web- 
type is shown by the red solid line in each panel. The mean and 
median cos(0) are also shown. If the two directions were randomly 
oriented with respect to each other the probability distribution 
would be uniform at unity. 95% confidence region of a random 
distribution is shown by the grey shaded region. All alignments 
are statistically very strong. 


filaments. V-web filaments and knots are clearly favoured 
by the Bisous model - these are over-represented (with re¬ 
spect to a random sampling of the V-web) by an extra 60% 
and 80% respectively. Sheets and voids are strongly unfa¬ 
vored, being seldomly selected by the Bisous model. In to¬ 
tal, 77% of Bisous filaments are either knots or filaments 
as defined by the V-web. These statistics indicate that the 
Bisous model and Shear tensor - although not in perfect 
agreement - agree well in terms of environment. It should 
also be added that this agreement occurs without any arbi¬ 



A1+A2+A3 = V' V ^ 6 

Figure 2. The sum of the shear tensor’s eigenvalues is the di¬ 
vergence of the velocity field (Ai + A 2 -I- A 3 = V ■ n ) which on 
the linear scales of the Wiener Filter is proportional to the over¬ 
density (V ■ V oc S). The distribution of the over-density in the 
four web types (voids, black dotted; sheets, green dashed; fila¬ 
ments, red solid; knots, blue dot-dashed) are shown here. Top 
panel shows this only for regions (cells) that are coincident with 
Bisous filaments, while the bottom panel shows this quantity for 
the full reconstructed CF2 volume. 

trary adjustment to the threshold used to segregate different 
V-web environments. 

The alignment between Bisous filaments and the 03 axes 
can thus be tested for environmental dependence by exam¬ 
ining how it changes as a function of web type. The four 
panels in Fig. 1 show the alignment for void, sheet, filament 
and knot regions. The mean and median cosine of the distri¬ 
bution is also indicated in each panel and is many standard 
deviations away from the value expected from a uniform 
distribution ((cos( 6 )) = 0.5). Although the alignment is sta¬ 
tistically very significant in all four V-web environments, the 
two methods agree best in voids. This is likely due to the fact 
that voids, being the most linear regions of the universe, are 
where the reconstructed velocity field most accurately rep¬ 
resents the true velocity and over-density field. 

The different environments (on scales of the Wiener Fil¬ 
ter) probed by our web classification is examined in Fig. 2 
where we show the distribution of the sum of the eigenval¬ 
ues (V • u = Al -b A 2 + A 3 ) for the Bisous filament locations 
tested in this study (top) and for the entire CF2 velocity field 
(bottom). The sum of the shear tensor’s eigenvalues is the 
divergence of the velocity which on the linear scales of the 
Wiener Filter is proportional to the over-density. Both top 
and bottom panels of Fig. 2 shows that knot regions (collaps- 
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ing along all three principal directions) are the highest over¬ 
density while void regions (expanding along all three prin¬ 
cipal directions) are the lowest over-density environments. 
If the distribution of matter in the universe constituted a 
perfect Gaussian random field, these four curves would be 
Gaussian distributions, whose standard deviation and mean 
increase as one goes from voids to knots. When the full CF2 
volume within 100 Mpc is examined (Fig. 2 bottom), the 
distribution of this quantity shows peaks in the sheet and fil¬ 
ament web types. Note that these two curves are fairly wide 
indicating that within the volume considered here, web types 
overlap in over-density (although the overlap is not extreme 
in the sense that for the most part overlap occurs only with 
adjacent web types.) A similar picture, but with different is 
seen in regions where the Bisous model is selected (Fig. 2 
top). However here the primacy of filaments and knots is 
clearly visible. In regions where the Bisous catalogue identi¬ 
fies filaments, these regions appear to be well sampled by the 
Wiener Filter reconstruction and are thus not by construc¬ 
tion the mean: as a result the eigenvalue sum distribution is 
considerably different than in the case of the full volume. 


4 CONCLUSIONS AND DISCUSSION 

There are many ways of looking at the large-scale distri¬ 
bution of matter. Traditionally this has been done by ex¬ 
amining the spatial distribution of galaxies obtained pri¬ 
marily through large redshift surveys. Such an approach, 
although intuitive, is susceptible to well known issues re¬ 
lated to redshift-space distortions such as fingers-of-god or 
the Kaiser effect. An alternative approach is to use the mea¬ 
sured velocity field combined with reconstruction algorithms 
to infer full (luminous-l-dark) underlying 3D density field. 
Although attractive because it probes the full 3D density 
field, this approach too suffers from well known problems of 
accuracy and precision due to sparse sampling. 

Nevertheless both methods can be used to study the 
so-called “cosmic web”, namely the complex network of fila¬ 
ments and voids that connect galaxies in the Universe. Here, 
we have studied two such methods: the “Bisous” model based 
on a marked point process and the “V-web” a dynamical 
classification of the large-scale structure based on the shear 
field. The former, results in a filament catalogue based on the 
2MRS galaxy catalogue while the latter uses a Wiener Fil¬ 
ter to reconstruct the velocity field from the Cosmicflows-2 
survey of peculiar velocities. Both methods have been used 
to examine the properties of galaxies within the cosmic web 
(see Gourtois et al. 2015; Metuki et al. 2015; Libeskind et al. 
2012, 2015; Tempel et al. 2015; Guo et al. 2015) 

The two methods, based on very different data sets and 
methods, agree well when it comes to identifying the prin¬ 
cipal directions of the large-scale structure. Fig. 1 shows 
that the spines of filaments from the Bisous model are well 


aligned with the axes that the V-web associates with fila¬ 
ments, es. This success is not trivial, because of the very 
different limitations of each individual method. A priori one 
would hope, but not expect, that these two methods agree. 

In fact, the agreement found here is considerably weaker 
than when the two methods are applied against each other in 
simulations, where biases and incompleteness play no role. 
Such a study was carried out by Tempel et al. (2014) who 
found a near perfect alignment between ea computed from 
the shear field and Bisous filaments (specifically, ~80% of fil¬ 
ament spines were aligned within 30° of each other). When 
applied to data, the number of alignments within 30° is 
closer to 23% - still well above the 13.4% one would expect 
for a random distribution, but considerably worse than in 
simulations. The reason the two methods fare poorer when 
applied to data is likely due to the Wiener Filtering recon¬ 
struction algorithm which smooths out the velocity field and 
perturbs its principal directions. Also, the disagreement may 
be a direct result of the fact that the two techniques are ap¬ 
plied to drastically different data sets. 

Perhaps the most important conclusion this study 
reaches, is a “proof of concept”. Reconstruction methods 
based on the velocity field appear to do a good job in recov¬ 
ering the filamentary network of the cosmic web, albeit with 
large effective smoothing kernels. In that respect, as mea¬ 
surements of the cosmic velocity field progress in parallel 
with deeper galaxy redshift surveys, the agreement between 
the two methods is also likely to improve. 
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